The dynamics of cracks in torn thin sheets 
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Motivated by recent experiments, we present a study of the dynamics of cracks in thin sheets. 
While the equations of elasticity for thin plates are well known, there remains the question of path 
selection for a propagating crack. We invoke a generalization of the principle of local symmetry to 
provide a criterion for path selection and demonstrate qualitative agreement with the experimental 
findings. The nature of the singularity at the crack tip is studied with and without the interference 
of nonlinear terms. 



I. INTRODUCTION 

This paper is motivated by some recent experiments in 
which thin sheets were torn under an out-of-plane shear 
mode (known also as mode III) [Tl|. The experiments 
tested a few initial crack configurations, including an ini- 
tial notch placed symmetrically in the middle of the thin 
sheet, but also richer configurations as shown in Fig. [T] 
and Fig. [51 The latter configuration is interesting in hav- 
ing two interacting cracks which appear to attract each 
other until a 'tongue' is separated from the mother sheet 
when the cracks coalesce. These experiments pose a chal- 
lenge for theory since it is not a-priori known what is the 
criterion for path selection for cracks in thin plates, or, 
in other words, where should the crack turn under the 
action of a given stress field. The aim of this paper is to 
address this question and to show that a proper combi- 
nation of elasticity theory in thin plates and a reasonable 
extension of the principle of local symmetry Q results in 
crack propagation that is in a reasonable agreement with 
the experimental results. 
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FIG. 1: A typical experiment discussed in Ref. 
(a) Schematic of the sheet before tearing : a centered cut was 
initiated in the film with inclined incisions at 45°. (b) Scan 
of the film after tearing : the path appears to be stable in an 
asymmetric configuration. 

The structure of this paper is as follows: in Sect. [ll] 
we remind the reader the well known theory of regular 
mode III fracture in a thick plate, and the modifications 



FIG. 2: Another typical experiment discussed in Ref. [H: 
Two cracks propagating simultaneously. The middle part of 
the film forms a tongue-like shape and detaches itself. 



required for thin plates. In Sect. Illll we discuss the nature 
of singular stress field at the tip of crack in a thin plate. 
Sect IIVI considers the modifications introduced by the 
large deflections that arc typical to the experiment under 
study. In Sect. |V]we describe our numerical procedures 
and compare the singular solution of the stress field to 
the theoretical expectations, including the influence of 
nonlinear terms. Next we consider the criterion of path 
selection and the resulting crack trajectories for various 
initial conditions. Comparison with the experiment are 
offered. Lastly, in Sect. I VI I we summarize the paper and 
offer conclusions. 



II. MODE III IN PLATES 

A. Thick Plates 

Imagine a crack which is situated along a symmetry 
plane (in the middle of a thick plate) subject to anti- 
symmetric distributed loads, directed in parallel with the 
crack edge. For a thick plate, we assume that the elastic 
displacements in the body are 



u = 0,v ~ 0, 
w — w{x,y). 



(1) 



The displacements in the xy plane are denoted by u,v, 
respectively, and the displacement in the z direction is 
denoted as w. Here and throughout this work we use the 
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Lagrangian convention in which displacements, strains 
and stresses are computed in the equihbrium position of 
a material point. According to Hooke's law, the compo- 
nents of the stress tensor are |3| 



— ^yy ^zz ^xy 
dw dw 



dy 



dx 



(2) 



where fj, = E/2{1 + v) is the shear modulus, E Young's 
modulus and v Poisson's ratio. 

In a cylindrical coordinate system r, Lp, z, where r = 
coincide with the crack edge and = ±7r with the crack 
faces, the stresses become 



dw 



1 dw 
r oip 



The condition of equilibrium is 



dr 







(3) 



(4) 



This equation can be written in a compact and 
coordinate-invariant form of the Laplace equation 



Aw = 0. 



(5) 



The Laplace equation is solved by any harmonic func- 
tions, with a specific solution selected by the boundary 
conditions. The crack faces are assumed to be traction 
free. 



dw 

dip 



= for p> = ±7r. 



(6) 



Due to the predominant anti-symmetry the harmonic so- 
lution that satisfied the boundary condition (jS]) can be 
expressed as the Fourier series Q], 



n=l,3,.. 



riLp 



(7) 



B. The Kirchhoff Thin Plate Theory 

In thin plates the result Eq. 1^] is no longer true. In 
order to analyze the correct stress singularity near the 
crack edge, we need to obtain the equation for the de- 
flection w (displacement in the z-axis) for a thin plate 
[^Q. We consider a plate of uniform thickness, equal to 
/i, and take the xy plane at the middle plane of the plate 
where 2 = 0. For a small deflection the displacements 
M, V vanish in this plane. However, the bending causes a 
small deformation along the z-axis, so the displacements 
become 



dw 
dx 



dw 
' dy 



(10) 



These create bending moments and a twisting moment 
that act on the plate elements and can be expressed as 



-D 



D 



f d'^w d^w 



\ dx"^ 

/ d'^w 
\dy^ 



dy"^ 

d'^w 
dx^ 



M. 



xy 



GxyZdz = D{1 — ly) 



dxdy 



(11) 



where D — Eh'^ /12{1 — i/^) is called the flexural rigidity 
of a plate. 

Substituting these expressions in the equation of equi- 
librium, the equation for the deflection w for pure bend- 
ing becomes 



A^w = 



(12) 



In addition to the bending moments and Aly and 
the twisting moment M^y, there are vertical shearing 
forces acting on the sides of the element. 



where An and Bn are constants. Applying this to Eq. 
yields the general form of the stresses 



n^l,3,.. 



(a„r"/2-i+6„r-"/2-i; 



/n(p\ 



n=l,3.,. 



cos ^-J (8) 



Barenblatt and Cherepanov |3| showed that the stress 
a^pz for small r has a singularity of the type r~^/^, and 
close to the crack tip it can be represented as 



Km 



IS — 

V2/ 



(9) 



where Km is the stress intensity factor which expresses 
the amplitude of the singularity. 



Qx 



(JyzdZ. 



(13) 



These forces can be expressed in cylindrical coordi- 
nates as I™ 



Qr = -D^{Aw), 
dr 



Qt = -D 



d{Aw) 
rdip 



(14) 
(15) 



For a crack in mode III, the crack propagation is con- 
trolled by the transverse shear stress a^z [1, 0| • Since it 
varies along the z axis, it is more common to take the 
sum of it along the thickness of the plate Eg. (ITS)) . In the 
next section we obtain the shear force Qt, and calculate 
the dominant terms close to the crack edge. 
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III. SINGULARITY OF THE SHEAR FORCE IN 
THIN PLATES. 

Imagine a crack cutting a thin plate along a symmetry 
plane, subject to anti-symmetric load. The crack tip is 
at r = and the crack faces are at (p = ±7r. For small 
bending, the deflection w is a solution of the bi-harmonic 
equation Ea. (|12l) . Assuming anti-symmetry, w can be 
expressed as the Fourier series: 



w 



^fe„(r) sm{nip) 



(16) 



where &„ (r) is defined by the boundary condition and the 
initial load. 

The crack faces are assume to be entirely free. Thus, 
along this edge there are no bending nor twisting mo- 
ments and also there is no vertical shearing forces f^, 



= Mr^ = = /or = ±7r (17) 

Kirchhoff Q proved that two boundary conditions are 
sufficient for the determination of the deflection w. Thus 
two on the above conditions are added together; in cylin- 
drical coordinates the transformed conditions at (p = ±7r 
read 



= -D 



1 dw 
r dr 
dAL 



1 d'^w d^w 



Q^2 



= 



= -D 



dr 
1 aAui 



Q3 



W 



1 d^w \ 
r drdif J 



(18) 
(19) 
= 0. 



The first boundary condition is satisfied identically when 
the Fourier series of Eq. (ITBl) is substituted in Eq. ((TH]) 
at if = ±7r. Repeating this with the second boundary 
condition (IT^ results in 



1 

D- 

r 



ldb„ 



^2 " jy_ " " 



Qj.2 
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ncos(±7m) = 
(20) 

This equation becomes zero at ip = ±tt when n ~ l/2+m 
and rTieZ. Thus, w can be expressed as 



2^ 6„(r)sin(-(^) 

Tl=l,3,.. 



(21) 



substituting this in Eq. p2)) . and the general solution of 
the deflection w takes the form of the following series: 



n=l,3 



2 ^-fL)„r2+t)sin(-^) 
(22) 

where An, BmCn and D„ are constants. 

The stress near the crack tip will be defined by the 
shear force Qt which relates to the anti-plane shear stress 
Tipz- From equations (|15l) and (1^^ . it follows that the 



shear force can generally be represented by an expansion 
of the form 



= J2 -^[a.(4-2n)r- 



Dn{4: + 2n)ri 



n=l,3,. 



n ,n . 



2 '2 

(23) 

notice that the first two terms of Eq. (1^^ give zero con- 
tribution to the shear force, since both of them are har- 
monic solutions for the Laplace equation. Now, taking 
the leading terms for n = 1, and we get that close to the 
crack tip the two dominant terms are 



= -Dicir- 



dir 



cos(f). 



(24) 



This result shows, that while in the standard mode III 
the inverse square-root is the dominant term close to the 
crack tip, for thin plates the dominant term scales as r~ = 
9] . We stress that this highly singular term is the result 
of linear theory. The numerical results presented below 
show that upon increasing the bending of the plate, this 
linear singularity is weakened due to non-linear effects. 



IV. LARGE DEFLECTION OF THIN PLATES. 

The theory of thin plates is applicable only for small 
deflections. In most cases it is applicable when the deflec- 
tion w is small compared to the thickness of the plate ^] . 
Thus, to relate to the experiments of Ref. jl| we must 
consider also large deflection effects. While for small de- 
flections the basic assumption was that the middle sur- 
face remains undeformed, here its deformation can not 
be neglected. Due to the large deflection, the quadratic 
terms in the strain tensor can not be omitted. For clarity 
of presentation we use now the notation u for the in-plane 
displacement with Ua and up replacing u and v of the pre- 
vious sections; we leave w to denote the z-component of 
the displacement. The strain tensor becomes 



dUa 



dXa 



1 dw dw 

2 dxa dxjs 



(25) 



The quadratic terms of Ua is neglected by assuming 
small deformation, but the same cannot be done with 
the derivative of w. Taking the nonlinear term into ac- 
count, the final equation for the deflection becomes 



t^a2 , 9 f dw 
DA w - h- — (Tap- — 

OXjS \ dXa 



0. 



(26) 



where a^p is the in plane stress acting in the middle plane 
of the plate. Assuming that there is no body forces or 
tangential forces acting in those direction we obtain the 
following equations of equilibrium: 



dxp 



0. 



(27) 
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The three equations ([26]) and ^ (Eq. ^ is actu- 
ally two equations), can be reduced to two equations by 
introducing the Airy stress function x defined by 



5^ 

9y2 ' 



dxdy' 



(28) 



Substituting in Eqs. (1261) and (|27|) we obtain the com- 
plete system of equations for large deflection of a thin 
plate 



^ I ox^ oy^ \ dxdy J 



2 ^ X d^w 
dxdy dxdy 



= 



= (29) 
(30) 



These equations are known as Fdppl von-Kdrmdn equa- 
tions. These equation are non linear and can not be 
solved exactly even in very simple cases. 

It should be noticed that the bending term (the first 
term of Eq. (|29|) ) has a higher power of h than the stretch- 
ing term {D ~ h^). Thus, for a thin plate subject to large 
deformation, the bending term can be neglected in com- 
parison to the stretching term and Eq. (|26p becomes 



her a/3 



dxadxp 







(31) 



In case of isotropic stretching, the stress components be- 
come the same in all directions, the equation for the de- 
flection becomes [1| 



Aw = 



(32) 



which is similar to the case of mode III for thick plates. 
The upshot of the last remarks is that for small deflection 
of a thin plate in mode III we expect that the anti-plane 
shear stress will result in a tip-singularity of the form 
r~2 . However, upon increasing the deflection, and there- 
fore increasing the in-plane stresses near the crack edge, 
this highly singular term should weaken, turning eventu- 
ally to a standard r~2 singularity. We will see that our 
numerical results support this expectation. 





FIG. 3: Color online. Left panel: the bi-Laplace solution 
for small deflection on thin plates under mode III conditions. 
Right side: The solution of the FvK equation which include 
the non-linear terms. The colors indicate the size of the de- 
flection w. 
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FIG. 4: Color online. A log-log plot of the shear force as a 
function of the distance from the crack. The three sets of 
data are for three different crack lengths. The best linear fits 
indicate the theoretical expectation that the stress close to 
the crack tip scales like r~ 
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V. NUMERICAL RESULTS 

A. The crack tip singularity. 

In order to analyze the stress field near the crack tip, 
we consider a system consisting of a rectangle thin elastic 
plate [L X L X h; h/L ^ 0.05]. An initial cut of length 
I was created at the middle of the upper edge. One side 
of the upper edge of the plate {—L/2 < x < 0;y = L/2} 
is bent to one direction and the other side {0 < a: < 
L/2; y = L/2} to the opposite direction, cf. Fig. |3l The 
other edges, including the crack's faces, are considered 
as free (cf. Eqs. [T7)) . We developed a solver using the 
Galerkin finite elements method on a triangular grid to 



calculate numerically the stress field and the shape of 
the plate. For the solution of the non-linear equations 
■ (PH)) ■ we resort to a combined incremental-iterative 
(Newton-Raphson) procedure. The solution of the shear 
force as a function of the distance from the tip for the 
bi-Laplace equation (fT2|) for a small deflection is shown 
in Fig. m We obtain that the singularity exponent of the 
shear stress close to the tip for several initial crack length 
is close to the theoretically expected value of -3/2. Next 
we investigate what happens upon increasing the bending 
of the thin plate. Close to the crack tip, the Gaussian 
curvature (the second term of Eq. (j30p ) is non zero and 
therefore an in-plane stress begins to accumulate. The 
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bi-laplace Eq. p2)) does not represent the situation any 
more; We must turn to the FvK Eqs. (HH) , ([201) ■ The 
difference between the solutions of these equations can 
be clear to the eye, as we demonstrate in Fig. [31 In the 
left panel the solution of the KirchhofF plate theory does 
not involve the creation of additional in-plane stresses. 
In the right panel the nonlinear theory creates in-plane 
stresses at the crack tip, resulting in a sizeable distortion 
of the plate shape ahead of the crack. This visual effect is 
also seen in the weakening of the stress singularity at the 
crack tip. In Fig. [5]we show, for several crack length, the 
singularity of the shear force. Note that for large bending 
the singularity of the shear force is weakened. 
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FIG. 5: The index of the shear force singularity r" as a func- 
tion of the amount of bending at mode III for several crack 
length. Increasing the bending creates an in-plane stress and 
weakens the singularity of the shear force at the crack tip. 



symmetrical about the crack direction. For thin plates we 
can also examine the subdominant term that breaks the 
symmetry for the deflection w. This term is proportional 

to cos{ip), 



cos{n(p). 
(34) 

substitute this in Eq. (jl5l) . and the first two terms vanish, 
being solutions of the Laplace equation. For n = 1 the 
third term also vanishes, and the first subdominant term 
that breaks the symmetry for thin plates comes from the 
last, term. Substituted in Eq. (fT5l) . it becomes identical 
to the term found in the Laplacian case Eq.(|33| 



Qt = -Dicir 



dir 



cos( — ) + Ai^2 sin((^). (35) 



We thus propose that the 'principle of local symmetry' 
in the present case dictates propagating the crack such 
as to annul A2. The generally accepted principle of lo- 
cal symmetry for crack propagation in mode I and II 
thus has a general analog for mode III as the condition 
A2 = Q. Of course, whether the crack propagates at all 
is dictated by the analog of the Griffith criterion, i.e. the 
energy release rate into the crack tip region from the elas- 
tic field must be larger than the dissipation involved in 
the creation of free surfaces. Note that in the standard 
case of 1/2 singularity one possesses analytic estimates 
of the energy release rate and by assuming a constant 
energy dissipation per unit area of free surface one can 
estimate the critical length of a crack. In the present case 
the analytic estimates of the coefficients ci and di in Eq. 
(|35|) are not available; their calculation requires further 
work, solving the bi-Laplace equation analytically. This 
calculation is beyond the scope of the present paper, but 
is planned for a future publication 



B. Crack Propagation 

The solution of the elastic problem, linear or nonlin- 
ear, does not contain a recipe how to advance the crack 
in a given stress field. For a precisely anti-symmetrical 
mode of tearing, the crack has no reason not to develop 
in a precisely straight trajectory. However, breaking the 
symmetry causes the crack to propagate in a different 
direction, and this direction is not known a-priori even 
when we know the stress field. We propose that in the 
present situation one can follow the strategy of Ref. |;3] 
which for the Laplacian problem introduced a subdomi- 
nant term to Eq.®. This term breaks the symmetry 

T^^ = '^^li. cos - sin((^) (33) 

The hypothesis of Ref. [3] was that a crack under lon- 
gitudinal shear propagates along a direction that annuls 
A2, making the stress distribution close to the crack tip 



C. Single crack trajectories with the principle of 
local symmetry 

We consider a plate of size 1x1 {—0.5 < x < 
0.5,-0.5 < y < 0.5}, with a small cut created in the 
upper edge {y = 0.5) at different locations parallel to the 
boundary. One side of the upper edge is bent to one di- 
rection and the other to the opposite direction (Fig [3|) . 
We chose the Young modulus E — 50 and the Poisson 
ratio ly = 0.3. We calculate the stress field for two cases: 
one for the Laplacian field, the normal mode III, and the 
second, for a thin plate, Eqs ((29 |) . (|30|) .To decide whether 
the crack should evolve or not we must choose a critical 
value of ci in Eq. p5p above which the crack must de- 
velop. To do so we computed the value of ci where the 
defiection equals the length of the initial crack, and chose 
this value of ci to be the critical value cl- This proce- 
dure is equivalent to choosing the energy dissipation per 
unit area of free surface when the Griffith criterion is em- 
ployed. In subsequent crack development we maintained 
the critical value of ci invariant. Thus the crack evolves 
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(a)Laplacian field (b)FvK solution 

FIG. 6: Color online. Crack paths for the normal mode III 
(upper panel) and for the thin plate in mode III (lower panel) . 
The thick line shows the initial crack. Note when the crack 
started in the middle of the upper edge, a; = 0, it continues 
straight for both cases, up to numerical inaccuracies. 




(a) Inclined cracks (b)Thin sheet with incisions 

FIG. 7: Color online. Upper panel: a comparison of the crack 
path of an inclined crack that starts at the center of the upper 
edge with an angle other than a right angle, and a crack that 
started with a right angle but not at the center of the edge. 
Lower panel: the calculated crack path with incisions of 45 
degrees inside the plate (red parallelograms). 



according to criterion ci > c* and the principle of lo- 
cal symmetry. Whenever the value of ci was below c*, 
we increased the bending quasi-statically by bringing the 
deflection closer to the crack's tip. Not surprisingly, we 
observed that a crack that was initiated in the middle 
of the upper edge a; = always continues straight. Fig. 
[SI However, when the initiation breaks the symmetry by 
starting the crack right to the center, the crack tends to 
turn to the right boundary. This tendency is stronger in 
the Laplacian case, where the crack curves sharply to- 
ward the right boundary. For the thin plates, the crack 
shows more stability by creating a smaller deviation from 
the straight line, and thus reaching the bottom in the 
present geometry. 

Next we investigated the trajectories followed by an in- 
clined crack. We initiate these with an angle with respect 
to the upper edge, such that a 'straight' initial crack has 
a right angle. In Fig. 7(a) we show that such initial 



cracks turn sharply towards the same trajectory that is 
followed by a 'straight' crack whose tip ends at the same 
point. This result indicates a surprising lack of memory 
for past history. The principle of local symmetry forces 
the inclined initial crack to turn sharply such as to agree 
almost exactly with the future trajectory of a straight 
initial crack whose tip is at the same point. Last, we 
compare the numerical solution to the experimental con- 
figuration Fig. [TJ We mimicked a similar set of incisions, 
and found that indeed the crack straightens up between 
the incisions as reported in the experiment. Even when 
creating incisions inside the sheet at an angle of 45°, the 
crack turns and continue almost in a straight direction 
(Fig. [7(b)). 



D. Interacting Cracks 




FIG. 8: Color online: the solution for the deflection w of the 
FvK equations for two cracks in mode III. The color code 
measures the deflection w. Notice that towards its end the 
sheet curves as a result of the coupling of the bending and 
the in-plane stress in the FvK equation. 
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The experiment showed in Ref. [l| showed that while 
a single crack is relatively stable, when two cracks propa- 
gate simultaneously in parallel to each other, they curve 
toward the center to create a tongue-like shape. In this 
simulation, we initiate two cracks symmetrically around 
the central axis, and bent the center to one direction, and 



FIG. 9: Color online. Numerical results for two cracks tra- 
jectories. Thick straight lines are the initial cracks, the dots 
are the results of the simulation; the thin line represent an 
exponent fit to the crack shape. In the inset, the numerical 
factor in the exponential function of the distance between the 
cracks (normalized to the width of the sheet). 
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the left and the right side to the other direction, Fig |8l 
The cracks evolve according to the conditions described 
in the last section. We allow both cracks to propagate at 
each step. By changing the distance between the cracks, 
we found that there exist a distance between the cracks 
at which both continue almost straight without attract- 
ing each other, cf. Fig [HI Theoretically one could expect 
that this minimal distance is 1/3 of the sample width, 
where the attraction to the edges cancels exactly the at- 
traction between the cracks. When the distance between 
the cracks is smaller, they attract each other. We found 
a fair analytic fit to the path of the cracks in the form of 
an exponential function y = a(l — exp{bx)), where x is 
distance from the intersection point of the two cracks, y 
is the distance between the two cracks, and a, h are con- 
stant. In the inset of Fig. IH] we present the numerical 
factor b ; when the distance between the crack is larger 
from l/L > 0.225 we can not find an intersection point 
and the numerical factor b is close to zero. 



VI. CONCLUSIONS 

In summary, we offered a criterion for path selection 
for cracks in thin plates. The criterion appears to gener- 



ate cracks that agree with those seen in the experiment. 
It therefore appears that numerical solutions can pro- 
vide reliable predictive tools for the dynamics of cracks 
in think plates in the Mode III tearing mode. In addi- 
tion we analyzed the stress field for mode III, for thick 
and thin plates, and showed that the singularity of the 
stress field close to the crack tip is stronger in the latter 
case, but it weakens when small in-plane deformations 
commence. 
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